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Abstract 

We propose a simple model where dark matter (DM) carries top flavor and couples to the 
Standard Model through the top quark within a framework of minimal flavor violation (MFV). 
Top-flavored DM can explain the anomalous top forward-backward asymmetry observed at the 
Tevatron, while remaining consistent with other top observables at colliders. By virtue of its large 
coupling to top, DM acquires a sizable loop coupling to the Z boson, and the relic density is set 
by annihilation through the Z. We also discuss contraints from current direct detection searches, 
emphasizing the role of spin-dependent searches to probe this scenario. 
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I. INTRODUCTION 

The nature of dark matter (DM) remains mysterious. Aside from its gravitational influ- 
ence, the particle physics properties of DM are largely unknown, and therefore it is worth- 
while to explore different avenues for how DM may couple to the Standard Model (SM). One 
interesting possibility is that the dominant couplings of DM to the SM arise through the top 
quark. As the heaviest known elementary particle, the top quark plays a fundamental role in 
many extensions of the SM, and experimental studies of the top quark are vital to validate 
our understanding of the weak scale. Anomalies in the top sector could be harbingers of 
physics beyond the SM, and this new physics may be connected to DM. 

At the Tevatron, both CDF and DO collaborations have measured an anomalously large 
top forward-backward asymmetry (Aps)- The measured parton- level inclusive Afb m the 
tt rest frame is (after unfolding and background subtraction) [TJ [2] 

J 0.164 ±0.045 CDF (9.41b- 1 ) 
FB ~ \ 0.196 ±0.065 D0(5.4fb- 1 ) 

for semileptonic tt events, in excess over the SM value Ap^ = 0.088 ± 0.006 [3] by ~ 2a. 
Notably, the discrepancy is larger for higher tt invariant mass (Mjj). For the high invariant 
mass bin (M« > 450 GeV), CDF has found Apf = 0.295 ± 0.066 [lj, a ~ 2.5a deviation 
from the SM value A^^' SM = 0.129^o;QQg. CDF has also observed a sizable Afb in dileptonic 
tt events, albeit with larger errors [4] . However, other top observables at both the Tevatron 
and LHC appear consistent with the SM thus far. 

In this work, we propose a simple model where DM, denoted Xt, carries top flavor and 
couples directly to the top quark. We take Xt to be the lightest component of a multiplet x 
carrying SM quark flavor, and we couple x to quarks according to the principle of minimal 
flavor violation (MFV) [5-8J through a scalar mediator <fr slightly heavier than t. DM stability 
is ensured automatically by MFV with no additional assumptions [9]. Although our model 
is phenomenological, we note that one realization of top-flavored DM is provided by warped 
grand unified models, where DM is a sterile neutrino partner of the top ^0TTl2] . (We note 
other related works exploring top-DM interactions [T3HT6] and flavored DM [TTl - fl9~] .) 

Top-flavored DM (TFDM) provides a number of interesting implications for collider phe- 
nomenology. Our model can explain Aps, while being consistent with other top observables. 
To generate Apb, an 0(1) forward-backward asymmetry is generated in uu — » <f)<p*, which 
is converted into an asymmetry in tt through decays — > txt- This mechanism was pro- 
posed for a supersymmetry- inspired model in ref. [20], but was later shown not to provide 
a viable DM candidate [21]. However, TFDM is a viable candidate for DM, evading many 
constraints associated with the model in ref. [20] as a natural consequence of MFV. Ad- 
ditionally, we study the collider implications of TFDM for the LHC, in particular the top 
charge asymmetry (Ac) and top-jet resonance searches [22]. 

On the cosmological side, TFDM can explain the observed DM relic abundance. As we 
show, DM acquires a one-loop coupling to the Z boson that is enhanced by the top mass. DM 
freeze out in the early Universe is governed by annihilation through the Z to SM fermions, 
providing the correct relic density for the same parameter region that can account for top 
Afb- TFDM in the local halo can give observable signals in direct detection experiments. 
In principle, the spin- independent (SI) DM-nucleon coupling can be large enough to explain 
the positive signals by CoGeNT [23] and CRESST-II [24], however there is strong tension 
with constraints from XENON [251 [26] and CDMS [23 M\ experiments. While SI signals 
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may be highly suppressed (through a Majorana splitting), TFDM predicts an irreducible 
limit on the spin-dependent (SD) cross section that cannot be evaded without additional 
channels to set the relic density. 

TFDM provides a minimal model to confront possible new physics in the top sector and 
DM simultaneously. Moreover, these different considerations independently point toward the 
same mass range for DM: m Xt ~ 20 — 90 GeV. In the remainder of this paper, we discuss 
each of these topics in detail. In Sec. [TTJ we present the TFDM model. We study top-related 



collider signatures at the Tevatron and LHC (in particular Afb) in Sec. Ill Relic density, 
direct detection, and other constraints are discussed in Sec. [IVj Finally, we conclude in 
Sec. El 



II. MODEL 



The principle of MFV ensures that new physics entering at the weak scale preserves 
the successful SM predictions of precision flavor measurements, which are sensitive to much 
higher mass scales [29]. The quark sector of the SM possesses a global flavor symmetry under 
G q = SU(3)q l x SU(3) UR x SU(3)d R that is broken only by Yukawa interactions. 1 Under 
MFV, it is assumed that new physics also preserves the symmetry G g , up to insertions 
of the Yukawa matrices Y u ,d acting as spurion fields. Although in most DM models, the 
DM particle is a singlet under G q , it is an interesting and less-explored possibility that DM 
transforms nontrivially under G q . One appealing feature of flavored DM is that DM stability 
is imposed automatically by the flavor symmetry [9] (for most but not all representations 
of G q ), whereas other DM models typically require additional symmetry assumptions for 
stability. 

In our model, we introduce two new fields, a Dirac fermion x an d a complex scalar 0: 

X ~ (1, 1, 0) SM x (1, 3, l) Gq , <f> ~ (3, 1, 2/3) SM x (1, 1, l) Gq , (2) 

where the two sets of numbers indicate quantum numbers under the SM gauge group 
SU(3)c x SU{2) L x U{1)y and G q , respectively. That is, x is a gauge-singlet, flavor multi- 
plet comprised of three states x — (Xu, Xc Xt), while is a flavor- singlet, color-triplet scalar 
with electric charge +2/3 which mediates interactions between x an d quarks. The Yukawa 
spurions transform as Y u ~ (3, 3, 1) and Y d ~ (3, 1, 3) under G q . 

The spectrum for the x states is constrained by MFV. Keeping only the top Yukawa yt, 
the mass term can be written as an expansion in powers (Y^Y u ) n with coefficients m n : 

- J? mass = x(m + mxY^Yu + ... )% = m Xu XuXu + m Xc XcXc + m Xt XtXt , (3) 

where the ellipsis denotes higher-order powers. Thus, Xu an d Xc are (approximately) de- 
generate, with mass determined by the zeroth-order term (m Xu ~ m Xc ~ mo), while the 
top-flavored state Xt has a different mass m Xt = (itlq + m\y\ + ...) from Yukawa insertions. 
Assuming these corrections are large and negative, we take Xt to be the lightest state, with 
mass m Xt rn Xu . 



We omit the lepton flavor symmetry SU(3)e L x SU(3) eR , which, though unimportant in our model, is 
relevant for lepton- flavored DM [TTUB] . 
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The scalar mediator couples x to up-type quarks qn = (u, c, t)n. As above, these 
couplings can be expressed in powers {Y^Y u ) n with coefficients g n : 

■2U = Qr(9o + QiYlYu + ■■■ )x<t> + h.c. = g u u R Xu<P + QcCrXc^ + 9ttRXi<t> + h - c - ( 4 ) 

Thus, according to MFV, g u ~ g c ~ go is given by the leading term, while g t = (go+giyf + ...) 
can be different due to higher-order insertions. We take these couplings to be real. 

The states Xu and 0, as well as the DM state Xt, play an important role in the phe- 
nomenology of our model. To summarize, the main parameters are the x masses (m Xu ,m Xt ) 
and couplings (g u ,9t) to quarks, as well as the scalar mediator mass m^. 

Next, we consider a possible Maj orana mass term for xt- 

/\ <YY), 

~£?Majorana = —^XfXt + h.C. (5) 

This term violates G q and cannot be accommodated within MFV, and therefore Xt must 
be Dirac. Effectively, if we regard the Majorana mass Am as a spurion, its representation 
is Am ~ (1,6, 1), which cannot be composed of Y u ,d- Nevertheless, it is useful to consider 
that small deviations from MFV may arise, and with Am ^ 0, the single Dirac state Xt is 
split into two Majorana components Xi,2, with mass m 12 = m Xt ± Am. In this case, even 
a small Majorana term Am <C m Xt becomes important for direct detection phenomenology, 
discussed in Sec. |IV| Clearly, we assume that deviations from MFV do not destabilize DM; 
this depends on underlying model-building assumptions that are beyond the scope of our 
phenomenologically- motivated study. 

Lastly, we discuss how the model of ref. [20] differs from TFDM. To explain Afb, ref. [20] 
introduces a scalar top ( "stop" ) partner t (playing the role of 0) and a single light neutralino 
X° (playing the role of Xu and Xt)- ^-fb arises by generating a forward-backward asymmetry 
in tt*, which is converted to ti through decays. However, having a single light state x° 
coupled to both -u^t and i^t allows for two decay channels t — > tx° and t — » ux°] the latter 
channel is phase space enhanced and dilutes A FB unless g t S> g u ~ 1- (Also, the charm 
interaction must be very suppressed due to D°-D° mixing bounds, g c j g u < 0.06 [2U]). This 
model is strongly constrained by LHC searches for jets plus missing energy (E^ ms ) [301 El] 
via pair production tt* — > uux°X°i and monojets plus E^ 1SS [32], [33] via ug — > tx° —> ux°X°- 
Also, identifying x° as DM is excluded by direct detection bounds, due to the large tree-level 
X°-u coupling [21]. For TFDM, all these issues are easily avoided by having x as a flavor 
multiplet. For m Xu > m^, decays 4> — > ux u and (ft — > cx c are kinematically blocked and do 
not dilute A FB since BR((fi — >• txt) = 1 without nonperturbatively large coupling g t . Signals 
from (mono)jets plus E^ 1SS are eliminated at tree-level since <fr ux u is forbidden (although 
loop-induced monojet signals can arise [IS])- Same-sign top production is not allowed by 
the Dirac nature of Xu, and no fine-tuning in the charm sector is required. Direct detection 



constraints are weakened since the DM xt couples to u at one-loop (see Sec. IV). 



III. TOP PHENOMENOLOGY AT COLLIDERS 

A. Top forward-backward asymmetry 

In the SM, the ti forward-backward asymmetry arises at order in the cross section, 
and therefore new physics in the top sector can provide a significant enhancement to Afb- 
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FIG. 1: (a): Aps is generated via production through a t-channel mediator (b): The dominant 
contribution to a negative Ac 



Several models have been proposed to explain the top asymmetry (e.g. axigluons 
and flavor-changing scalar or vector bosons [I2H5H]), mostly relying on interference between 
new physics (NP) and gluon amplitudes in ti production. However, the simplest models are 
strongly constrained by top observables at the Tevatron and LHC [4*0 1 1601461] . as well as by 
atomic parity violation constraints [65], requiring additional model-building ingredients to 
be viable. 

For TFDM, Afb is generated entirely through NP without interference with the QCD 



la 



amplitude, shown in Fig. 1(a), through the mechanism of ref. [20]. First, an 0(1) forward- 
backward asymmetry is generated in 00* through t-channel Xu exchange, and then this 
asymmetry is transmitted to tt through decays — > txt and 0* — > ixt- Fig- 
parton-level differential cross section 

da(uu 00*) _ fit giPl{l - cos 2 6) 



gives a 



dco$6 32tts (l + /35-2^cos0 + 4m2 /s) 2 



(6) 



where (3^ = y 1 — 4m 2 i /s, which is strongly peaked in the forward direction (cos# > 0) as 

long as 4m 2 u /s < 1. Also, since QCD production of scalars 00* is p-wave suppressed, this 
channel does not give a large contribution to inclusive tt production. Searches for tt plus 
E™ 1SS (see e.g. [SSHSH]) are evaded by assuming to be nearly degenerate with txt (i-e., 
m<l> — rn t — m Xt w 0). We also take m Xu > to maximize Aps by having BR(0 — > txt) — 1 
(since — > ux u , c-Xc are forbidden). 

To explain the Apb data, the TFDM model must generate a large inclusive asymmetry, 
as well as an asymmetry that rises with M t j, without introducing large corrections to the 
total inclusive ti cross section (a t t) or the differential cross section (da / dM t A . The inclusive 
ti production cross section is a t i = 7.65 ± 0.42 pb [691471] . consistent with the theoretical 
SM predictions of cr| M = 7.067 ± 0.26 pb [72], allowing for cx^- p = (0.58 ± 0.49) pb at la. 

For da/dM t t, we quantify this constraint in terms of the ti cross section for M t i > 450 
GeV, measured by CDF to be aJJ gh = 1.92 ± 0.48 pb [73], while the SM prediction is 
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Mgh,SM _ 2. 17 ±0.10 pb [71]. This observable, as we show, provides the strongest constaint 
on our model. We have also investigated the high M t i tail, quantified in terms of the second- 
highest invariant mass bin (Mti = 700 — 800 GeV) following refs. |1HI [7S]- However, within 
our model, this observable is not as constraining as 0"|j gh . 

To simulate this signal, we generate NP events (pp — > tixtXt) at the partonic level in 
MadGraph5 vl.4.3 [76J, with a pp initial state at yfs = 1.96 TeV and m t = 173.3 GeV [77J. 
We compute the total (SM + NP) A FB as 



fT S - M a N - p 
A _ A SM °tt , 4 NP °tt f 7 \ 

FB — rt-FB _SM+NP ~+~ FB _SM+NP ' V ) 



S_M+NP 1 "-FB SJVI+NP 
a tt a tt 



where A F ^ is the asymmetry generated in ti produced from NP-only. We work to leading- 
order (LO) in computing NP cross sections; we take SM quantities given at their highest- 
available order. We take A F ^ to be the most recent next-to-LO (NLO) QCD prediction 
including LO electroweak contributions j3j I78TI80] . We expect that the inclusion of (un- 
known) NLO corrections to NP quantities would likely reduce the allowed parameter space, 
and therefore our treatment is conservative in the sense of allowing for the largest possible 
parametric region. 



B. Top charge asymmetry 

At the LHC, a related observable is the top charge asymmetry (Ac) with respect to the 
boost direction of ti pair. The ATLAS and CMS collaborations have measured inclusive 
asymmetries (at 7 TeV) [8TL 182] 

f -0.018 ±0.036 ATLAS (1.04 fb- 1 ) , , 

c ~ \ 0.004 ±0.015 CMS (5.0 fb- 1 ) ' W 

consistent with the SM prediction, A C M = 0.0123±0.0005 PES]. Additionally, the combined 
inclusive ti cross section at the LHC is a t i = 173.3 ± 10.1 pb [83], compared to the NNLL 
SM calculation, cx| M = 162.6 ± 17.1 pb [M] for m t = 173.3 GeV. 

There are two contributions to the charge asymmetry in our model, shown in Fig. [TJ and 
they have opposite sign. The positive contribution comes from the process uu — > (fxfi* — > 
tixtXt- The negative contribution arises from ug — > which in turn decay to Xu uixt 
and (f) t~xt- The incoming w-quark is harder than the incoming g at the LHC. As a result, 
the outgoing i is harder than the outgoing t, resulting in a negative Ac- The cross section 
for the CP conjugate process is an order of magnitude smaller, due to the incoming u parton 
distribution function. As a result, our model can generate a positive Afb at the Tevatron 
while also generating Ac of either sign at the LHC. In particular, the process mediated 
by Xu(p generates a sizable Ac ~ —20%. This effect, first discussed in [851 ]86j , can be 



understood from the parton-level cross section computed from Fig. |l(b) 

daiug ->• (pXu) gig 2 



dcos9 32ns 24 



(l-/^cos0), (9) 



neglecting terms of order {m 2 ^ — m 2 ^)/s. This process tends to produce anti-aligned with 
the initial u, giving a negative contribution to Ac- 
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To compute the charge asymmetry in our model, we generate the above two NP processes 
at the partonic level in MadGraph5 vl.4.3, with a pp initial state at y/s = 7 TeV. Similar 
to ApBi we compute 

A, , - ASM °tt | ANP °tt (in\ 

X-C - SM+NP +/1 C SM+NP' ^ iU J 

°tt °tt 

where is the asymmetry generated in ti produced from NP-only. 



C. Top+jet resonance 

The heavy \u can be searched for in top+jet resonance searches at colliders [221 E3 EE]- 
This signal arises in tij events in t-channel mediator models; the mediator is produced on- 
shell, in association with a t, and has a two-body decay to ij. Within our model, a top+jet 
resonance arises from sequential 2-body decays Xu u 4> — > uixt, where Xt is undetected. 
The top+jet invariant mass is given by 



m tj — m t H ^2m^ + m t ~ m xt ~ cos @\J ( m l ~ m t ~ m X t) 2 ~ 4m t m x t j > I 11 ) 

where 9 is the angle between the i and u in the (f) rest frame. Since (f) decays isotropically, the 
top+jet resonance has a "box" feature. However, since collider constraints prefer a squeezed 
spectrum (such that m Xu — m Xu and — m t — m xt ~ 0), the width of the box 

is typically smaller than the jet energy resolution [SH]. In this case, we can approximate 
mtj ~ m Xu — m Xt , which is valid within 5% over the parameter space of interest. 

The ATLAS search [HE] places the strongest limits on the production of top+jet res- 
onances within the context of specific models. For the color-singlet vector model {W), 
the parton- level cross section is constrained to be less than (23, 14, 7) pb for mw' = 
(200,300,400) GeV, respectively. Assuming the acceptances and efficiencies are compa- 
rable, we apply these limits to TFDM, with the replacement mw — > m Xu — m Xt , and we 
extrapolate between these different mass points following [88J. We do not extrapolate below 
200 GeV, since searches [53 EH] have not been performed in this range. 



D. Analysis 

In this section, we identify what is the region of parameter space where TFDM can 
account for an anomalous Apg, while remaining consistent with other observables. We 
simulate the NP contribution to top observables in MadGraph. We emphasize that the top 
asymmetry generated in TFDM does not necessitate a light Xt-, unlike the case in [20J. By 
restricting to a degenerate spectrum, m^ — m t — m Xt w 0, we can increase the mass of both 
Xt and 4> while minimizing EJp lss in ti production. 

Our main results are presented in Fig. [2} In general, the collider phenomenology of our 
model depends on five parameters (g u , g t , m Xu , m Xt , m^), and the different panels in Fig. [2] 
illustrate different slices through this parameter space. In panel (a), we show our results in 
the (m Xt ,g u ) plane, assuming = m Xt + 175 GeV and m Xu = + 10 GeV. In panels 
(b— d), we show our results in the (m Xu ,g u ) plane, for specific choices of m xt . We assume 



7 




m Xu (GeV) m Xu (GeV) 

(c) (d) 

FIG. 2: Allowed la region that fits top observables at Tevatron and LHC in the (m Xt ,g u ) plane (a), 
and the (m Xu ,g u ) plane with m Xt = 45, 60, 80 GeV, respectively (b,c,d). Inclusive Aps (A^ h ) 
is consistent with the blue (red) region, with the overlap of both shown in purple. Remaining 
Tevatron (LHC) observables are shown by the solid (dashed) lines, with allowed regions indicated 
by arrows. The hatched region is consistent with all observables at la except c!f gh - 
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9t = 9u for simplicity; since gt affects only the total width of 0, this choice is of minor 



importance. 2 

The different top constraints are shown as follows: 

• The shaded red band is consistent with A^^ = (29.5 ± 6.6)%, while the shaded blue 
region is consistent with inclusive Afb = (17.4 ± 3.7)%, taking an error-weighted 
average of Eq. ([I]). The central shaded purple band shows the overlap consistent with 
both. To be conservative, we consider that Ap^ and A^!| h,SM are at the upper limit of 
their la allowed values. 

• The solid lines indicate Tevatron constraints. The solid green line shows the la region 
for the inclusive ti cross section aJf Y . The red contours indicate a^l gh,]<lp at the level of 
0.49 and 0.74 pb, corresponding to a 1.5a and 2a enhancement over the difference of 
observed and SM values (aJJ gh,NP = -0.25 ± 0.49 pb). The arrows denote the allowed 
region. 

• The dashed lines indicate LHC constraints at 7 TeV on the inclusive ti cross section 
a t| HC (green), top+jet resonance searches (orange), and Ac = (0.1 ± 1.4)% (brown), 
taking an error- weighted average of Eq. ([8]). The arrows denote the allowed region. 



The hatched region is consistent with all aforementioned measurements at la except for 
tension with cr.j g . Taken at face value, the preferred (hatched) region for TFDM exhibits 
some tension at the level of 1.5 — 2a with cr^ lgh . However, reconstruction and acceptance 
efficiencies, as well as higher-order corrections, which we have not addressed, may be impor- 
tant. 



Fig. 2(a) shows the allowed DM mass and coupling that is consistent with Tevatron and 
LHC top data. A DM mass m Xt ~ 20 — 90 GeV is allowed, with coupling g u ~ 1.4 — 2.2. 
Tevatron ti measurements provide the strongest constraints, while m Xt > 90 GeV is excluded 
by Ac and m Xt < 20 GeV is excluded by cr^ HC . It is noteworthy that relic density constraints 
also point to a similar DM mass range, discussed in Sec. |IV| 

Figs. 2(b) , 2(c) , 2(d) show how the allowed parameter region depends on m Xu , for m Xt = 



45,60,80 GeV, respectively. At larger m Xu , the parameter region becomes more restricted 
by Tevatron ti measurements. Improvements in sensitivity in top+jet resonance searches 
may play an important role in exploring TFDM parameter space. 



IV. RELIC DENSITY, DIRECT DETECTION, AND OTHER CONSTRAINTS 

Top-flavored DM, by virtue of its large coupling to top, naturally acquires a significant 
one-loop coupling to the Z boson. This coupling sets the DM relic density and also generates 
observable signals in direct detection experiments. From the diagrams shown in Fig. |3j the 
resulting effective interaction is 

2L t = a z —xa lt PLXtZ lt , (12) 
cw 



In Sec. IV we show that relic abundance considerations fix g t as a function of 
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FIG. 3: Feynman diagrams for DM coupling to the Z boson. 



where gi is the SU(2)l gauge coupling, cw = cos 6^ is the cosine of the weak mixing angle, 
and az is the new physics coefficient, given by 

Sg? mf 2 2 . . x — 1— log a; . . 

az = 32^^ K/m ^' FZ{X) = (1-x) 2 ■ (13) 



Although Eq. (12) corresponds to a dimension-six operator Xtl^ 1 PiXtH^ D , the mass scale 
suppression 1/iV is compensated by the Higgs vacuum expectation value (H) 2 . Within 
Fig. [3j (H) manifests as the mass of the top quark in the loop. If the new physics scale 
is comparable to mt, as we have assumed, there is no mass scale suppression of the DM 
coupling the Z, giving az ~ O.Olgf . 

In the early Universe, DM annihilates to SM fermions // through the Z boson. The 
total annihilation cross section is 

working at leading order (s-wave) in the relative velocity v, except we keep the full -in- 



dependence of the propagator to properly treat resonant annihilation [HO]- In Eq- (14), 
s ~ 4m^ t (l + 1> 2 /4) is the center-of-mass energy, sw = sin^vy, and N*, Qf, and T 3 f denote 
the number of colors, electric charge, and weak isospin, respectively, summed over all SM 
fermions / excluding the top, which is not kinematically allowed. 

For a Dirac particle, the required value is (av) an k 6 x 10~ 26 cm 3 /s to achieve the ob- 
served DM relic abundance during freeze-out. 3 Although (crv) an is suppressed by the loop 
suppression of az, it is enhanced by resonant annihilation for m Xt ~ rriz/2, by the large 
number of SM fermion final states, and by a large coupling g t > 1. Following ref. [90J, the 
DM relic abundance is 

_ l2 1.03 x 10 9 GeV _1 

n dm h 2 « — - 2 — , (is) 

g*s/g* m p iJ(x fo ) 

where X{ = m Xt /T{ (Tf D is the freeze-out temperature) is given by the solution to the 
equation Xf Q = log(0.038g Xt m p i (crt>) an / y/g*Xf ), with g Xt = 2. Post-freeze-out annihilation is 
parametrized by the integral 

1 f°° dx I x 3 r°° 2 

J ( x fo) = o / -^Han, {<ru) m =\ — dvv 2 {av) an e- xv /4 , (16) 

2 Jx to % V 4vr J 



To clarify our notation, we note that (av)& n is the thermal average of (av) a 
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FIG. 4: Feynman diagrams contributing to the SI DM-nucleon cross section relevant for direct 
detection experiments. For each case we show only one representative diagram, while we include 
all possible diagrams in our analysis. 



where the 1/2 arises for Dirac fermions. We evaluate these integrals numerically, which is 
important on- resonance (m Xt ~ m^/2), where the s-wave approximation breaks down. 

Next, we discuss several experimental constraints. If DM is lighter than m^/2, decays 
Z — > XtXt are allowed, contributing to the invisible Z width. The partial width is 

T{z Xttt) = tScJ7 (1 ~ m2jml) v 1 - 4 ^'^ 1 w L6 MeV x (S) 2 • (17) 

The measured invisible Z width is consistent with the SM, Y{Z — > inv) expt — T(Z — > inv)sM = 
— 1.7±1.5 MeV [91j, and any new physics contribution must be no greater than O(MeV). We 
also note that TFDM does not contribute to highly-constraining parity- violating observables 
at one-loop, unlike standard t-channel models [22] . 

Direct detection experiments, which search for nuclear recoils from DM in the local halo, 
have important implications for top-flavored DM. In particular, the DAMA/LIBRA |92j, 
CoGeNT [23], and CRESST-II [23] experiments have found evidence for a signal from 
C(lOGeV) mass DM; however, taken at face value, these observations appear at odds with 
null results from the XENON100 [23], XENON10 j2g], and CDMS [27] EE] experiments. 
Reconciling these different results in an area of active debate, with appeals to astrophysical 
assumptions, experimental systematic effects, or DM model-dependent assumptions (see [93T 
196] and refs. therein). 

Although Xt has no tree-level couplings to nucleons, couplings arise at one-loop order, 
shown in Fig. |4j For a given nucleus N, the SI DM-nucleon scattering cross section is (see 
e.g. [93(98]) 



2 



a 



si /% 



N - — (Zf p + (A-Z)f n ) 2 (18) 

7T 

where (Z,A) are the proton and atomic numbers, and /i^v = m N m X t/ ( m N + m xt) * s the 
reduced mass. The proton/neutron scattering amplitude is f PtTl = + / p b ° x + f£ n + , 
where the different terms delineate each contribution in Fig. [4} The Z and box terms provide 
the leading contributions to <Tat through an effective vector-vector operator qj^qXtJ^Xt, with 
q = u,d. The Z vertex contribution is 

ti = q jr< /p z = -a-^)^f> (is) 

and the box contribution is 

/r = ^f, fr = a -f, a b o X = I |^^o X KK,<K) (20) 



ii 



with loop function 



. _ x 2 (y- l) 2 logx -y 2 (x - l) 2 logy+ (x - l)(y - l)(x-y) ( . 

^ boxMj - 4(x-l) 2 (y-l) 2 (x-y) ' 1 J 



The Higgs-exchange amplitude is 



eh = G F gjm v , n m Xt 
' p ' n 167r 2 V2m 2 



* , (22) 



where « 0.83 is a hadronic matrix element [HH]- For the gluon amplitude, using results 



' TG 

from ref. |lUUj . we obtain (for the contribution denoted therein as fa 



fln = -%rn xt m p , n f& n \r + + ? + ) , (23) 



I 
36' 

with loop functions f+ l defined in ref. [100] (with quark mass m t and "squark" mass m^). 
The Higgs-exchange and gluon amplitudes are quantitatively much smaller than the Z and 
box terms, providing 0(0.1%) corrections to f p>n . The 7 exchange contribution corresponds 
to a magnetic dipole interaction, with a Xt dipole moment 

eg 2 m x 2 2 1 + 2x log x - x 2 
fel 7(m(/ ^' F 7 (x) = • (24) 

For 10 GeV DM, a magnetic dipole moment at the level of 10~ 18 ecm can potentially al- 
leviate tensions in direct detection results due to the momentum dependence of the cross 
section [101J. Here, we have fi xt ~ 7 x 10~ 21 g 2 ecm for pa m t , which is too small unless 
g t in nonperturbatively large, in which case f% is similarly enhanced. Therefore, we neglect 
H Xt in our analysis. 

In Fig. [5j we show how constraints from direct detection experiments map onto top- 
flavored DM parameter function of (g t ,m Xt ). For definiteness, we have fixed 

— m t — m Xt = 2 GeV; we have also taken g u = 2 and m Xu = 400 GeV for evaluating / p b ° x , 
although these parameter choices are of minor importance since the box term provides only 
a subdominant contribution to the scattering cross section. The solid black line (labeled 
Qdm) shows the parameters required for Qdmh 2 = 0.11 [91], while the gray region is excluded 
for T(Z — > XtXt) < 3 MeV. Provided Xt is a Dirac state, constraints on SI scattering provide 
the strongest limits, due to the sizable one-loop coupling to the Z. The left panel of Fig. [5] 
shows that the entire parameter space that gives the correct relic density is excluded by 
XENON100, as well as other measurements. These constraints require m Xt < 10 GeV or 
g t < 0.1. In this case, annihilation through the Z is insufficient to achieve the correct relic 
density, and additional channels would be required. 

The stringent limits from SI scattering can be evaded if DM is Majorana. We suppose 
that the Dirac state Xt acquires a flavor-breaking Majorana mass term, given in Eq. (j5j), 
that splits Xt into two Majorana states X12, with only the lighter state xi populated in the 



Universe today. In terms of X12, Eq. (12) becomes 



^int = clzt^- (X2l^Xi + Xil^X2 + Xi7 M 75Xi + X27 M 7sX2) Z„ . (25) 

The vector interaction becomes inelastic. If the mass splitting Am is larger than (9(100 keV), 
X\ has insufficient kinetic energy to access the excited state X2 |103j . On the other hand, 
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Spin-independent bounds for Dirac DM Spin-dependent bounds for Majorana DM 




2 5 10 20 50 100 2 5 10 20 50 



m Xl (GeV) m Xl (GeV) 

FIG. 5: Left: For Dirac DM, we show exclusion curves for SI scattering from XENON100 [25] 
(blue), XENON10 [26] (purple), and CDMS [271 EE] (brown dashed), while enclosed regions show 
signal regions from CoGeNT [23] (green) and CRESST-II [21] (red). Right: For Majorana DM, 
we show exclusion limit on SD scattering from XENON100 [102J, as well as the prospective reach 
for limits for SD scattering on neutrons (solid brown) and protons (dashed brown). Gray region is 
excluded by T(Z — > inv). Black line gives the DM relic density Qdmh 2 = 0.11. 



the axial vector interaction remains elastic, contributing to the SD DM-nucleon scattering 
cross section. Following ref. [97], the SD cross section on the neutron or proton is given by 



°5 = (d«Aj»»> + d d Af n) + d s A^) (26) 

for axial- vector interaction Jz? cff = dq^^TsSXiTjuTsXi- The coefficients d q are given by 



d d = d s = , (27) 



2V2 8 ' " A 2^2 



due to Z exchange and box contributions (d u only). We take A„ = Ajj = 0.78, 
Ai n) = -0.48, Ai p) = Ai n) = -0.15 [99]. 



Fig. [5] (right) shows how constraints from SD DM-nucleon scattering map onto TFDM 
parameter space. XENON100 provides the strongest limits on TFDM scattering on the 
neutron |104j . shown by the blue line (see also ref. [102J). 4 The combination of the constraints 
from XENON100 SD scattering and T(Z — > inv) exclude g t > 5, while a large parameter 
region remains viable. The brown solid (dashed) contours indicate the SD cross section (in 
cm 2 units) for DM scattering on the neutron (proton), showing how future searches will 
impact TFDM. The relic abundance constraint provides a lower bound > 10~ 43 cm 2 , 
which is saturated for m Xt ~ 45 GeV. Future experiments, such as XENON1T |105j . at the 
level of <t^ d ~ 10 -41 — 10 -42 cm 2 offer the potential to explore nearly all of TFDM parameter 
space, except for very near the Z pole. 



4 Constraints from SI interactions (via Higgs and gluon amplittudes) are much weaker and lie above the 
vertical scale of the plot. 
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The Majorana mass splitting Am can also play an important role for DM annihilation. 
The annihilation cross section for X1X1 ~ >* // is p-wave or chirality-suppressed, given by 

<"-> = 8^(( 3 -S + m |r|) S Nl { Ti < m > + l my(T >> ~ 2M + *»*>) • 

(28) 

while X1X2 annihilation has an unsuppressed s-wave component given by Eq. ( Jl4| ). This 
implies that DM annihilation is suppressed in the halo today (since only xi is populated), 
and thus indirect detection signals are quenched. However, if Am < Tf D ~ GeV, then 
X2 is populated during freeze-out, and the correct relic density is obtained by X1X2 f f 
coannihilation, as in the Dirac case. Thus, the relic density line (black) in Fig. [5] (right) 
provides a lower bound on g t . The bound is saturated for Am <C T fo , while larger Am ~ T fo 
requires larger values of g t since the X1X2 annihilation rate has a Boltzmann suppression 
exp(-Am/T fo ). 

Although we have thus far assumed a standard (symmetric) freeze-out cosmology, it is 
also possible that Xt might be asymmetric DM (see |106[ I107j and refs. therein). Within 
extra-dimensional realizations, TFDM carries a generalized baryon number, and therefore 
the dark sector may have a baryon asymmetry similar to visible sector (in fact, they may 
be naturally related) [TOldT]. To the extent that Am « and Xt~Xt annihilation is efficient, 
the relic density is set by the initial dark baryon asymmetry. The natural mass scale for 
Xt is (fidm/^b) x m P ~ 5 GeV. While there is no Majorana splitting to evade SI limits, 
asymmetric DM naturally favors the low mass window where SI limits are weakened. 5 On 
the other hand, efficient symmetric annihilation requires (ca>) an > 6 x 1CT 26 cm 3 /s. Since 
annihilation through the Z is insufficent in the low mass window, additional channels are 
required in this case. 

V. CONCLUSIONS 

The top forward-backward asymmetry remains a persistent anomaly from the Tevatron. 
Many new physics models have been proposed, with new degrees of freedom introduced 
mainly to address this anomaly. Here, we have shown that possible new physics in the top 
sector may be connected to another long-standing particle physics puzzle - the nature of 
dark matter. We have considered a simple phenomenogical model where dark matter Xt 
carries top flavor and couples to the SM via the top quark within the framework of MFV. In 
particular, we find that top physics and dark matter constraints independently point toward 
the same mass range for DM: m Xt ~ 20 — 90 GeV. 

For collider phenomenology, the scalar mediator plays a key role, with an O(l) forward- 
backward asymmetry in on-shell cjxp* production transferred into a ti asymmetry through 
decays — > txt, as in ref. [20]. Missing energy signatures in ti events are avoided for 
— m t — m Xt ps 0, and several constraints on the model of ref. [20] are evaded due to the 
additional flavor structure of our model. We have shown that top-flavored DM generates 
a sizable top asymmetry and is consistent with other Tevatron and LHC observables for 
m Xt ~ 20 — 90 GeV and for a perturbative range of couplings, albeit with some tension the 



5 Introducing a Majorana splitting as large as O(100 keV) leads to Xt~Xt oscillations that wash out any 
initial asymmetry, leading to the usual symmetric freeze-out scenario. 
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high invariant mass ti cross section measured by the Tevatron. Future top+jet resonance 
analyses should be able to explore the parameter space of the model. 

The DM relic density and direct detection signals are determined by an effective coupling 
of DM to the Z boson. This coupling arises at one-loop and is enhanced by the mass 
of the top appearing in the loop. The DM freeze-out abundance is set by annihilation 
XtXt ~^ ff through the Z vector coupling. We have shown that the correct relic density can 
be achieved for DM masses prefered by collider constraints (especially near the Z pole). If 
DM has a small (flavor-breaking) Majorana mass splitting 100 keV < Am < GeV, direct 
and indirect detection processes are governed by the Z axial-vector coupling. Stringent SI 
limits are evaded, and the most promising direct detection signals are for SD scattering 
(with a lower bound cr?° > 10~ 43 cm 2 at m Xt « m^/2). Annihilation in the halo today is 
p-wave or chirality suppressed. Our model illustrates that DM annihilation through the Z 
boson remains a viable mechanism for fixing the relic density, and future SD experiments 
are crucially important for exploring this possibility. 
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